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ON A GENERAL METHOD FOR TREATING TRANSMITTED MOTIONS 
AND ITS APPLICATION TO INDIRECT PERTURBATIONS* 

BY 

ERNEST WILLIAM BROWN 

The mathematical treatment of any physical problem demands the construc- 
tion of an ideal problem in which the conditions are different from those of the 
actual problem. It is assumed that the same general laws hold for the actual 
and ideal problems, but tlie complexity of the circumstances surrounding the 
former makes simplifications of some kind necessary in order that the analysis 
should not be unreasonably tedious. The ideal problem is therefore usually 
constructed by neglecting at first some of the influences which form a part of it 
but which are assumed to affect the results to a much smaller extent than those 
we retain. The simplified problem, which I call problem A , is then solved. 
The second step consists in finding what changes are necessai'y in the solution 
of problem A when some or all of the neglected influences are included ; this 
second problem I call B. The question under consideration here is the deduc- 
tion of the solution of B when that of A has been found. 

One of the methods for solving B is that known as the Variation of Arhi- 
trary Constants. In the deduction of the solution of A , arbitrary constants a 
arise, and numerical values, depending on the given conditions of the problem, 
can be assigned to them. In order to obtain the solution of B, it is assumed 
that the solution of A is still available provided that a different set of values be 
assigned to the constants a . But the new set of values will usually be different 
at different times and will therefore depend on the time ; in other words, the 
constants a become variables in problem B. 

This last fact gives us another view of the relation between A and B. Sup- 
pose that the position of the system which constitutes B is defined by variables 
X, and that we desire to change from the variables x to other variables 2^ in 
order to simplify the treatment of B . Then one p-set, namely, the a-set, can 
be obtained by supposing that the x-set and the ^-set are connected by the same 
equations which connect the x-set and the a-set in the solution of problem A . 

* Presented to the Society under a different title, February 25, 1905. Received for publica- 
tion April 8, 1905. 
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Thus A is no longer looked upon as a dynamical problem but as a collection of 
relations which give certain transformations useful for the treatment of B. 

But we can go further than this. There is no special necessity imposed on 
us to use any of the a-sets thus found. We may be able to discover another 
set of variables, say a ^)'-set, allied to the cc-set by i-elations different from 
those which connect the a;-set and the a-set — relations suggested indeed by the 
latter and perhaps possessing similar properties, but possessing also additional 
properties more useful for the purpose in hand. 

It is proposed to investigate here one such ^'-set. The relations between the 
X- and a-sets usually contain a certain number of other constants a which were 
present in the original equations of problem A. It is assumed that the x- and 
p -sets are connected by relations of the same form as those which connect the 
X- and a-sets, but that the a have given variable instead of given constant 
values assigned to them in the ^'-set. Certain advantages are gained in the 
application of the method to the consideration of the 'indirect inequalities' in 
the motion of the moon, i. e., in the consideration of those perturbations of the 
lunar motion which arise from the action of a planet on the earth's motion and 
which are transmitted to the moon on account of the predominating influence of 
the earth on the moon. I have given some indications of this ^'-set in a former 
paper.* Unfortunately, the main idea was obscured by the method adopted to 
develop it, and some statements concerning the applicability of it to the lunar 
problem require further limitations than those previously imposed. The present 
method of treatment also biings to light a theorem (art. 10 below) which I 
believe is new and which is of considerable importance in the discussion of the 
indirect inequalities. 

1. The method of Jacohi. The equations of a dynamical system are sup- 
posed to have been expressed in the canonical form 

dx, d<b dy. dd> 

where, to fix ideas, the 7i variables x^ may be thought as the coordinates defining 
the position of the system, the n variables y. as the corresponding momenta, and 
<f> as the sum of the kinetic and potential energies; these interpretations are 
nevertheless not necessary. 

Let the variables t x, y he changed to 2w new variables^?, q, connected with 
X, yhy the relations 



* On the variation of llie arbitrary and given constants in dynamical equations, TransaotiODS 
of the American Mathematical Society, vol. 4 (1903), pp. 333-350. 

fThe sufiSxes will be omitted whenever no confusion is thereby caused. In all cases t, j 
take the values 1 , 2, • ■ • , n . 
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dS dS 

(2) 2'^ = ax/ ^' = 5^;. \:s=s(^,p,t)-\, 

then Jacobi has shown that the equations satisfied by the new variables are 

in which, after dSjdt has been formed fi;om S expressed in its original form as 
a function of t and the x, p, the function <}> + dS/dt is supposed to be ex- 
pressed in terms of ^, §', « by means of (2). No assumption is made regard- 
ing the form of S, but in order to secure the independence of the x, y and also 
that of the i? , 5' , it is necessary that the determinant 



d'^S 



(i,j=l,2, •••, »), 



dx.dp. 
shall not be zero. 

Suppose that the x, y receive two sets of independent arbitrary variations 
Sx, hy and 8'x, S'y, inducing corresponding variations in the^J, c[. Equations 
( 2 ) give 

Whence, if the x , y are expressed in terms of t and the p, q, and if we adopt 
the usual notation : 



we obtain the following relations : 






(5) U'.'l,}-^^ 
We also have 


[i>o?J = 


(»+i), 


[?., ?J = o 


(i=or +i). 



d4> d 

^ = -4>(:c,y,t). 

2. Particular forms of 8. Suppose first that S is so chosen as to satisfy the 
equation ^ + dS/dt = . The^, q are then constants and they may be con- 
sidered as the arbitraries of the solution. Such forms of S are solutions, con- 
taining n arbitrary constants exclusive of that additive to S, of the partial dif- 
ferential equation 

Next suppose that S is so chosen as to satisfy the equation 
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dS 

(6) ^+al = ^' 

where 5 is a function of the p only ; this is a form much more convenient in 
many of the applications. Then the p and therefore the q are constants. Let 
us put 

(7) i',= c., q. = w.= b,t+e., h,= -~^-, 

i 

the arbitrary constants being c, e, and B being a function of the c only; at 
present B is undetermined but it may be determined so as to satisfy future con- 
ditions as to the forms oi x, y, when the latter are expressed in terms of c, e, t. 
I shall suppose that this problem has been solved and that this second form 
for S has been used so that the solution is expressed in terms of c^ , w^, t,a\l of 
which are independent. On account of (2), we have, 

dS dS 

(8) yi=dx' ^"i^W lS = S{x,e,t)-i, 

i i 

[C;,w.]=l, [c.,w.] = (i^j), 

(9) 

[c.,cj = 0, [w.,Wj] = (i=or+i). 

3. The constants present in <f>. Suppose that (f) contains m constants a^ to 
which numerical values have not been assigned. Some or all of these constants 
will be also present in every form oi S. I propose to investigate in this article 
certain properties of S and <f) with relation to these constants. It is not neces- 
sary to suppose that the a^ are all the constants present in (/>. Such of them 
as are not specified in the a, can be supposed to have had numerical values 
assigned to them, either explicit or implicit, so that they do not enter into the 
discussion. 

In order to make the following results more general, I suppose that the m 
constants a,^ have been replaced by tn variables or constants u^^ connected with 
the a^ by the m independent relations 

(10) /,(«,,«„ = {h,k = l,2,---,m),*- 

so that the u^ are functions of the a,^ and t only. 

The function S was originally a function of the x, c, a,t and it now becomes 
a function of the x, c , u, t; denote it by S' when it is expressed in the latter 
form. We have 

4 114 ft 



*The values 1,2, - • •, m are given to A, k in all cases. 
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Put 
Then 



(^■^^ dt ~ ■^' du, dx. ^' ^ ^* Bu, du, ""* ■*" du^ dt ' 

Also denoting by Q' the expression of any function Q in terms of the x , c 

u and t, and putting 

d'<p d 

/i h 

we obtain 



du^ - aw, * dy. du,^ - du,^ \dt '^ ) •""' du,8a. 

d / dS' dS' „\ „ . d\ 



, dx, 



du. dx. ' 

h I 



Combining with (11) we find 

*^^^'' dt ~ du ^""^"du, du, ' 

ft ft ft 

an equation which enables us to find l/,^ directly. 
Again 

d S' dS dxj 

'' ~ du, ~ du. ~ ^^' du, ' 

/( h It 

where S and the x are expressed in terms of the c , w, u, t. Hence if t/, be 
similarly expressed, 

dU,^ d^S d dXf 

^ ' dc. ~ dc.du. dc. ^^'du, 

1 t h t h 

But since the c, u, w are independent, 

d'S d'S 



dc.du. du, dc 



'. "" du, \ dc. i dx. dc. ) 

% ft \ I 3 t / 

d ( „ cZx,\ d dx, 

= d^X''+^^yrcrcJ = dii:^jyjTc: 

Substituting in (13), performing the differentiations with respect to m,, and 
remembering the notation (4), we find 
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(14) ^^=K,cJ- 

i 

In an exactly similar manner it can be shown that 



(15) ;77r=[«.-«^J> 



dw. 



~dM^~ du^ ~ L^''' "*-■■ 

The Uj^ can therefore be determined directly from the expressions for x, y 
when the solution has been completed. 

4. Two limitations. In the remaining portion of this paper, I suppose 

1°. That t is not present explicitly in ^{x, y, u, t), so that if t was present 
explicitly in the original expression for </> , one or more of the u must be variable. 

2°. That t is not present explicitly in the x, y when the latter have been 
expressed in terms of the c , w, u, and t . 

These assumptions do not very much limit the scope of the results and they 
render the exposition rather more simple. 

We have now 

(16) ^ = 0, ^^ = s,C7;a, 4>^\u,u, = -b. 

5. Variation of arbitrary constants. The usual method consists in sup- 
posing that, after the equations (1) have been solved in terms of the c , w , and t , 
a new function J3(x, y, t) is added to <^ ; account is taken of i? by supposing 
that the x, y retain the saxoB forms when expressed in terms of the c, w, and t, 
so that the c, e are no longer constants. This amounts to nothing more than a 
transformation from the old variables x, y to the new variables c, w, and the 
results of art. 1 can therefore be used. 

In fact we now have 

^ ( dx. dx, . dx\ d 

Adc.^dio.^dtJ dx.^^ " 

which, on account of (9), give 

dcf d Q dwj d Q 

dt^dw.' 'dt^~^/ 
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where Q = cf) -{■ H + dSjdt is supposed to be expressed by means of (8) in 
terms of the c, m), and t. But being a function of the x, y, and t., and 
d Sjdt a function of the x, c, and t, and neither nor d S/dt containing total 
derivatives with respect to t, 4> + d S/dt must retain the same form whatever 
be the values given to the c, w. Hence we still have (f) + dS/dt=— B,b, 
function of the c only, and 

Q=-B + B. 

Suppose that i? is due to changing one of the constants a present in (^ to a 
variable. Then 

d'6 d'^d) 

where Sa is the variable part added to a. The same method may be pursued, 
it is true, but if we follow it we take no account of the fact that the ybrm of (p 
considered as a function of the x, y, a, t is unchanged, and it would seem as if 
this unchanged form should be made use of to simplify, if possible, the results. 
The new method of varying the constants which follows recognizes this fact and 
brings to light certain results which will be of value in the applications. As in 
art. 3, the a will be replaced by the u ; the two conditions of art. 4 will also be 
imposed. 

6. Variation of arbitrary and given constants. I suppose that the x, y 
retain the same forms when the m in are replaced by functions of the 
time different from those with which the original problem is solved, so that we 
replace not only the c, w but also the u by their new values in order to get the 
new values oi x,y . The new values of c , w> will, therefore, be different in this 
method from those found by the method of art. 5. In fact we now have, since 
t is present in the x , y only through its presence in the c , w , u, 

_, ( dx. dx. \ _^ dx. d(f) 



On account of (9), (14), (15), these give 



'' dn. * dx. 



(17) 

/ J J 
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if Uj^ be independent of c, w., and 

The equations are therefore still canonical. 

Now when the u, c, w have their original values, denoted by the index , 

(18) </>» + 2,C^«<=-5». 

This equation is an identity after the x, y ^ ii, f7have been replaced by their 
values in terms of the w, c , u, and t , and we can, therefore, give to the c,w,u 
any arbitrary variations, symbolically denoted by 8^, so that 

(19) </," + \4> +\{Ul + B,U,){ul + 8,«J =-B'>-S^B. 

The actual variations in Q are produced by adding portions, denoted sym- 
bolically by S , to the c , 'w, u, so that 

(20) Q = </." + 8<^ + 2„(C/;. + SC/,)|(w« + 8«„). 

As <j>, U, B, u are functions of the c,vi,u only and not of their derivatives 
with respect to t, we can take h^ = h. But we cannot in general put 
S^u = dBu/dt. For the identity (19) demands that u be first expressed in 
terms of u by means of (10) and the variation then taken, while in ^ a portion 
Su is to be added to u and the derivative with respect to t then formed. It 
was the failure to recognize this distinction that caused the error in my former 
paper.* The method there used implies the substitution oi — B for Q under 
the false assumption that the derivative of Su with respect to t would be the 
same as Bii. The error is immediately evident in the case of those u^^ which 
were originally constants ; for them, m,_ = and, therefore, Sii^^ = , while obvi- 
ously dhujdt will only be zero if hu be constant, and this will not generally be 
the case. 

Combining (19), (20), with 8j = 8, we obtain 

(21) q=-B + ^,U,(J^^hu,-hu)j, 

the final form of Q . 

I have also obtained this result, when squares and higher powers of S are 
neglected, by direct transformation from the ordinary canonical equations of 
art. 5 above, making use of the relations given in art. 3. 

7. There is, of course, no special necessity to confine this method to those 
variations which are produced by varying the given constants in (^. Suppose 

*Loo. cit., art. 17. 
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that the present method is retained for the consideration of a function S added 
to (f), while the given constants in <^ do not vary. We must then add to Q the 
function 

the first term of which is necessary to counterbalance the variations of the u in 
(f> . But this proceeding appears to be of little or no value, as the results are 
more complicated than in the usual method. On the other hand the variations 
of the u are now entirely arbitrary and we might be able to so choose them that 
the form of Q is more simple or the canonical equations more easy to solve. 
What advantage may be gained will depend on the particular problem in hand 
and I shall not pursue this development further in the present paper. 

8. Case of the m,^ satisfying caiionical equations. Suppose that the u are 
divided into two classes c^ , w^ = h'jt -\- e^ , where i^ = — dB' jdc[. I suppose 
that in the initial problem, the c' , e' are constants and that i?' is a function of 
the c' only. Instead of using the symbols U, put 

C" = — — . W = [S'=S(x.,c.,c'.,w'.)-\, 

h h 

so that 

(22) ^ = - 5 + 2,. q 1 8c; + 2, tf; (| k - ss;) . 

Further, I suppose that the variations of the c' , w' are given by 

or 

d dR' d dR' 

dt '' d%!}\ ' dt " 5c„ * 

Substituting in (22), we obtain 

where it is supposed, in accordance with the previous work, that i?' is indepen- 
dent of the c, ^o. 

9. Properties of the solution in the case of the lunar theory. The applica- 
tion of these results to the lunar theory is immediate when the necessary limita- 
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tions which constitute this particular case are imposed. For the sake of brevity, 
I shall not prove these limitations here but merely state those of them that are 
necessary for my purpose. They are as follows : 

1°. The function R' is the disturbing function due to a planet on the motion 
of the earth. It can be expanded as a sum of periodic terms and it has a small 
factor whose square can and will be neglected. 

2°. A periodic term in R' of period a years will give rise to periodic terms 
in the 8c' of the same period with a factor a in the coefficients: in the 8w' the 
factor d' also occurs. 

3°. It is supposed that the orbit of the earth lies in a fixed plane. There are 
then two angles w^ , and three angles w^ . 

4°. The W'l^ consist only of sums of cosines of multiples of the four differ- 
ences of the five angles w,^, w^, the terms independent of t being zero; the C^ 
consist of sines with similar arguments. 

5°. The c\^ are functions of the mean motion and of the eccentricity of the 
earth's orbit. 

10. The indirect planetary inequalities in the nioon^s motion. There are 
very few sensible planetary inequalities in the moon's motion except when a is 
a number fairly large compared with unity ; only the cases where a is large are 
treated here. It is further supposed that squares of the small factor which 
accompanies R' are neglected. Then from (17), (23), 



where 



(K 5 (>, c? ^ so, ^ db. - ^ db. . , 



If R' contains a term of long period, a years, the factor a may be intro- 
duced in the 8c, 8c' and therefore the factor a^ in the 8w, 8zo', and no higher 
power of a can be present. Hence the general 

Theorem. When squares and higher powers of the ratio of the mass of a 
planet to that of the sun are neglected, the large factor a dtie to a long period 
inequality can never occur in the corresponding term in the moons motion to 
a power higher than a^, even if its square Is present in the corresponding 
inequality of the earth's motion,. 

The importance of this result lies in the fact that it gives us a means of 
rejecting in advance many terms of long period, which under the old method 
would have to be examined in order to find out whether their coefficients were 
sensible in the moon's motion. The advantage of this method lies, to a great 
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extent, in the theorem just stated and in the calculation of particular classes of 
inequalities ; this will appear in the following article. 

11. Special classes of inequalities. I divide the long-period inequalities into 
three classes for which certain results can be immediately deduced. The order 
of these refers to the power of a present. In the periodic inequalities it is 
supposed that at least one of the arguments belonging to the disturbing planet 
is present. 

(o) Inequalities containing the lunar arguments. They are of two kinds. 
First, those produced by terms of moderately long period in the moon-earth-sun 
system combined with terms of nearly the same period in the earth-sun-planet 
system. The former terms have all quite small coefficients and there will be 
few of such inequalities with sensible coefficients. Second, those produced by 
combining terms of short period in the two systems. We shall generally be 
able to neglect the he' , hw' , and the principal parts of the coefficients arise only 
in the Zw^ obtained from limiting the equations to 

dt ' dw.' dt ' ■' dc. J 
(6) Periodic inequalities not containing the lunar arguments. Here 

Sc. = 

7 

and the C , W are limited to those terms which have the arguments w^ only. 
The latter have all small coefficients and moreover the terms produced by Q^ 
have only the factor a. Hence it will be possible to obtain the chief parts of 
the coefficients by solving the equations 

d db. , 

, ^^i = 2. -rj Be. , c. = const. 



dt ' * dc 



(c) Non-periodic inequalities. These again are of two classes. First, the 
terms containing t to the first power only. These merely give constant addi- 
tions to the mean motions of the lunar arguments. Second, the secular inequal- 
ities, which contain t^. These latter are determined by 

8c. = , -rr Bw. = S, -,-; Be, , 

' dt ' ''del 

where Bc',^ is determined from the non-periodic term of H'. This is Newcomb's 
theorem.* The theorem therefore only applies exactly, with the given limita- 

* Aciion of the Planets on the 3Ioon, Washington Astronomical Papers, vol.5 (1894), 
p. 191. 



1905] AND INDIRECT PERTUKBATIONS 343 

tions, to the secular inequalities, but it is approximately true for the greater 
part of the inequalities independent of the lunar arguments, while it is not 
applicable to inequalities containing the lunar arguments. In arts. 17-19 of 
my former paper (loc. cit. ante), it was stated that the theorem applied to all 
ndirect inequalities. 

I have contented myself here with a brief outline of the chief results which 
the method furnishes for the lunar problem. There are numerous details and 
properties connected with the functions C , W which might find a place. But 
these are of more interest in connection with a numerical calculation of the 
inequalities than of theoretical interest in the variation of the constants or even 
in the problem of three bodies. It therefore appears advisable to give them 
later with the publication of the details of a method for such an investigation. 

Haveefoed College, 
April 8, 1905. 



